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The two-dimensional attractive Hubbard model is studied in the weak to intermediate coupling 
regime by employing a non-perturbative approach. It is first shown that this approach is in quanti- 
tative agreement with Monte Carlo calculations for both single-particle and two-particle quantities. 
Both the density of states and the single-particle spectral weight show a pseudogap at the Fermi 
energy below some characteristic temperature T* , also in good agreement with quantum Monte 
£NJ ' Carlo calculations. The pseudogap is caused by critical pairing fluctuations in the low-temperature 

renormalized classical regime {hui < ksT) of the two-dimensional system. With increasing tempera- 
ture the spectral weight fills in the pseudogap instead of closing it and the pseudogap appears earlier 
D ' in the density of states than in the spectral function. Small temperature changes around T* can 

5^ , modify the spectral weight over frequency scales much larger than temperature. Several qualitative 

■ results for the s— wave case should remain true for d-wave superconductors. 

71.10.Fd, 71.27. +a, 71.10.-w,71.10.Pm. 
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O ' I. INTRODUCTION 

,o,; 

For the past several years pseudogap phenomena found in the underdoped high temperature superconductors |l| 
and organic superconductors || have attracted considerable attention among condensed matter physicists. For these 
materials the low frequency spectral weight begins to be strongly suppressed below some characteristic temperature T* 
that is higher than the transition temperature T c . In the high temperature superconductors, this anomalous behavior 
has been observed through various experimental probes such as photoemission specific heat ||, tunneling ||, 

NMR |Q, and optical conductivity Although various theoretical scenarios have been proposed, there is no consensus 
at present. These proposals include spinon pair formation without Bose-Einstein condensation of holons j^-|ll[], stripes 
|^2~|l5|, hidden d— density wave order strong superconducting fluctuations [p7|-p5|, amplitude fluctuations with 
dimensional crossover J2(| and magnetic scenarios near the antiferromagnetic instability ]27|]28| ] . 

Although the above theories for the origin of the pseudogap are very different in detail, those that do not rest on 
spatial inhomogeneitics can be divided, roughly speaking, in two broad categories : Weak-coupling and strong-coupling 
explanations. In the strong-coupling approaches, the single-particle spectral weight is shifted to high-energies. There 
is no weight at zero frequency at half-filling. That weight however increases as one dopes away from half-filling, as 
qualitatively expected from the Physics of a doped Mott insulator. Recent angle-resolved photoemission spectroscopy 
(ARPES) experiments (2^,0 find, in the superconducting state, a quasiparticle behavior consistent with this point 
O ■ of view. If we consider instead a weak-coupling approach, either in the strict sense or as an effective model for 
quasiparticles, the only known way of obtaining a pseudogap is through coupling to renormalized classical fluctuations 
in two dimensions. In the repulsive two-dimensional Hubbard model in the weak to intermediate-coupling regime, 
analytical arguments J3l],[52| and detailed Monte Carlo simulations Q strongly suggest that indeed antiferromagnetic 
fluctuations can create a pseudogap in the renormalized classical regime of fluctuations. This mechanism has been 
confirmed recently by another approach [B4| but earlier studies had not found this effect. j35|,|36| 

The present paper focuses on superconducting fluctuations and the attractive Hubbard model in weak to intermedi- 
ate coupling. The purpose of the paper is twofold. First, in section [n], we validate, through comparisons with Monte 
Carlo simulations, a non-perturbative many-body approach [B7J that is an extension of previous work on the repulsive 



model [ |3l| , |32| . Formal aspects of this method are presented in the accompanying paper. |37| Then, in the second part 
of the present paper (section [n|) we study the mechanism for pseudogap formation due to superconducting fluctua- 
tions. ]3q| More extensive references on pseudogap formation in the attractive Hubbard model may be found in section 



III. Earlier Monte Carlo work I39J37I] and analytical arguments f33,p9[ have suggested the appearance of a pseudogap 



in the renormalized classical regime of pairing fluctuations. We study the appearance of the pseudogap in both the 



1 



density of states and the single-particle spectral weight A( k p,ui), showing that, in general, they occur at different 
temperatures. General comments on the relation to pseudogap phenomena in high-temperature superconductors may 
be found in the concluding paragraphs. 



II. A NON-PERTURBATIVE MANY-BODY APPROACH COMPARED WITH MONTE-CARLO 

RESULTS 

In the first subsection, we present our approach in simple terms. More formal arguments are in an accompanying 
paper p7[ . In the second subsection, we show that our approach is in quantitative agreement with Monte Carlo 
simulations for both single-particle and two-particle quantities. 



A. A non-perturbative sum-rule approach 

We consider the attractive Hubbard model for electrons on a two dimensional square lattice 

where = — 2t(cosk x + cosk y ), U is the on-site attractive interaction (U < 0) and N the number of lattice sites. 
Throughout the calculations, the constants t, H, ks and lattice spacing are taken to be unity. The index a represents 
spin. This Hamiltonian is not a valid model for d-wave superconductors, but it is the simplest model for which 
it is possible to check the accuracy of approximate many-body results against Monte Carlo simulations. Once the 
accuracy of the many-body technique has been established, it can be generalized to the d— wave case. Furthermore, 
many qualitative results do not depend on whether one has s— wave or d— wave pairing. 

The non-perturbative approach to the attractive Hubbard model presented in the accompanying paper is an ex- 
tension of the approach used in the repulsive case |3^]. In the first step (which was called zeroth order step in the 
repulsive model case) , the self-energy is obtained by a Hartree-Fock-type factorization of the four-point function with 
the additional constraint that the factorization is exact when all space-time coordinates coincide. It is important 
to note that this additional constraint, analogous to the local field approximation of Singwi et al. leads 
to a degree of consistency between one- and two-particle quantities that is absent from the standard Hartree-Fock 
factorization. Functional differentiation, as in the Baym-Kadanoff approach |jt2| , then leads to a momentum- and 
frequency-independent particle-particle irreducible vertex that satisfies |43| 



(1 -n T )(n x )' 



U PP = U )Y ^ . (2) 



With this approximation, the particle-particle susceptibility, which obeys the Bethc-Salpeter equation illustrated on 
the first line of Fig.[j], can now be written as, 

where the irreducible particle-particle susceptibility is defined as 

Xo%) = §E G ^- fc ) G -^( fc )- (4) 

k 

The vertices and the Green functions in xp 1 ^ (?) ar e at the same level of approximation in the sense that the irreducible 
vertex U pp is obtained from the functional derivative of the self-energy entering The vertex U pp = (5T,/8G) is a 

constant and EW in zero external field is also a constant, leading to a Green function that has the same functional form 
as the non- interacting Green-function G°(fc), where k = (^ik n ,}tj stands for both the fermionic Matsubara frequency 

ik n and the wave vector ft. The constant self-energy E^ 1 ) can be absorbed in the chemical potential by working at 
constant filling. If needed, we have argued in the accompanying paper |3^] that the following should provide a useful 
approximation for £W since it allows the first-moment sum rule for the pair susceptibility to be satisfied |37|], 
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FIG. 1. The first line is a skeleton diagram representation of the Bethe-Salpeter equation for the susceptibility in the 
particle-particle channel (Eq.(46) of accomp any ing paper |$^) and the second line is the corresponding equation for the 
self-energy (Eq.(58) of accompanying paper |37|). In the Hubbard model, the Fock contribution is absent, but in general it 
should be there. Solid lines are Green's functions and dashed lines represent the contact interaction U. The box and attached 
lines are the particle-particle susceptibility while the ellipse is the irreducible particle-particle vertex. The particle-particle 
susceptibility is obtained by identifying points 1 and 3 in the Bethe-Salpeter equation. 



s (i) _ U _ U pp C 1 Z g) (5) 

At this first level of approximation, only the double occupancy (D = (n-j-nj.)) is needed to obtain the irreducible 
vertex U pp and hence the pairing fluctuations. The value of (ri|ri|) can be borrowed from some exact calculations 
or approximate estimates but, as in the repulsive case, we found that accurate results are obtained when (n^n^) 
is determined self-consistently from the following "local-pair sum rule" , a consequence of the fluctuation-dissipation 
theorem for the s-wave pairing susceptibility 

§ Xp(«) expH<z„Cr) = (At A) = (n T n 4 ) . (6) 

Substituting Xp ] Eq.(|) for the pair susceptibility and Eq.(|J) for the irreducible vertex U pp leads to an equation that 
determines double-occupancy, and hence U pp , self-consistently 



N i i tt ((i-"tKi) v (i)/^ 
? L + u (i-„ T >( ni >Xo (1) 



exp(-i<7„0 ) = (n^n^ . (7) 



This first part of the calculation is referred to as the Two-Particle Self-Consistent (TPSC) approach. |44| 

Once the pair susceptibility has been found, as above, the next step of the approach consists in improving the approx- 
imation for the single-particle self-energy by starting from an exact expression where the high-frequency Hartree-Fock 
behavior is singled out. This expression is represented by skeleton diagrams on the second line of Fig.(|l|). The piece of 
the exact self-energy that is added to the Hartree-Fock behavior represents low-frequency corrections. These involve 
Green functions, vertices and pair susceptibility for which we already have good approximations from the first step 
of the calculation. One thus substitutes in the exact expression the irreducible low frequency vertex U pp as well as all 

the other quantities at the same level of approximation, namely (k + q) and Xp^ (o) computed above, obtaining 
in Fourier space 

^(k) = Un^ -U^J2 U pP 4 1) (q)G W a (q - k), (8) 

1 

where q = (iq n ,lf) stands for both the bosonic Matsubara frequency and the wave vector. Here T is the absolute 

(2) 

temperature. The resulting self-energy T, a (k) on the left hand-side is at the next level of approximation so it 
differs from the self-energy entering the right-hand side. Physically, it is a self-energy coming from taking into 
account cooperons. As stressed previously |3^|||, it is important that the irreducible vertex U pp , (or T^) as well as 
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G [ a '(k + q) and Xp (<?) all be at the same level of approximation otherwise some results, in particular with regards 
to the pseudogap, may come out qualitatively wrong. 

The particle-particle irreducible vertex U pp may be regarded as the renormalized interaction strength containing 
vertex corrections. Note that in the expression for the self-energy E^ 2 \ Eq.(||), one of the vertices is bare while the 
other one is dressed with the same E3j particle-particle irreducible vertex function U pp that appears in the paring 
susceptibility. In other words, we do not assume that a Migdal theorem applies. If this had been the case, the self- 
energy would have, amongst other things, been proportional to U 2 instead of UU PP . In T-matrix theory the bare U 
is used everywhere which, in particular, leads to a finite-temperature phase transition at the mean-field temperature. 
In our case, as described below, we proceed differently, avoiding altogether a finite-temperature phase transition in 
two-dimensions. 

We briefly summarize some of the constraints satisfied by the above non-perturbative approach |3j| . In Eq. (Q) , A is 
the local s— wave order parameter Ci[C^. Anti-commutation relations, or equivalently the Pauli principle, imply that 
( [A, A^] } = 1 — n. This in turn means that the convergence factor exp(— iq n 0~) in the local-pair sum rule is necessary 
because ( [A, A^] ) = 1 — n implies that, except at n = 1 , one needs to specify if r = + or r = 0~ in the imaginary-time 
pair susceptibility. Either one of these limits however leads to the same value of U pp since our approach |37j satisfies 
exactly this consequence of the Pauli principle [jl6): ([A, A^] ) = 1 — n. In complete analogy with the repulsive case 
discussed in Ref. [[IT]], one can invoke the 2D phase space factor, 2nqdq, and the Ornstein-Zernicke form of the pairing 
correlation function near a critical point to show the following. Deep in the renormalized classical regime, where 
the characteristic frequency v c of the retarded pairing susceptibility satisfies v c < T, the superconducting correlation 
length £ increases exponentially £ ~ exp(C /T) with decreasing temperature. In the latter expression, C may be 
temperature dependent. This behavior is the one expected in the O (n — oo) universality class |^7|] and not in the 

XY, or O (2) , universality class, where £bkt ~ exp(C / (T — Tbkt) 1 ^ 2 )- The precise dependence on temperature of 
the correlation length in the temperature range between the beginning of the renormalized-classical regime and the 
actual critical regime is not known analytically. In the regime that we explore, the pseudogap begins to open at a 
temperature T that is quite a bit larger than Tbkt- The latter was estimated to be at most of order of O.li for \U\ = 4 
by Moreo et al. pq |. As in the repulsive case, we will see below that the pseudogap appears in A( k p, uj) if the pairing 
correlation length grows faster with decreasing temperature than the single-particle thermal de Broglie wavelength 
ith = vf/T. 



B. Comparisons with Monte Carlo calculations 

In this section we show, by comparing with Monte Carlo calculations, that the present non-perturbative ap- 
proach is an accurate approximation. The calculations are performed for the same lattice size as the correspond- 
ing Monte Carlo calculations. It is important to note that at half-filling the Lieb-Mattis canonical transformation 

Cil —> exp (— iC$-~f*ij cj|, with Q = (7r,7r), maps the attractive model onto the repulsive one, pair fluctuations at 

wave vector if being mapped onto transverse spin fluctuations at wave vector if + T$ . With the proviso that in 
the attractive model at half-filling one would need, because of symmetr y [p3| ]5lfl , to take into account the charge 
fluctuations, we can state that the comparisons done in the repulsive case ]32|] |33|] apply for the canonically equivalent 
attractive case. We note in particular that it was shown that the convergence to the infinite-size limit is similar in 
the Monte Carlo and in the non-perturbative approach p3[ | . We restrict our discussion to cases away from half- filling. 
The Quantum Monte Carlo simulations [j52[ |p3fl that we performed were done using a Trotter decomposition with 
increment At ~ 1/10 in imaginary time and the determinantal approach ]54j |. Typically, about 10 5 or more Monte 
Carlo sweeps of the space-time lattice are performed. 

We begin with double occupancy. Table I shows the value of double occupancy calculated from Monte Carlo 
simulations for various temperatures and fillings. The last column shows the value obtained from the self-consistent 
equation for double-occupancy Eq.(Q). Clearly the accuracy of the approximation deteriorates at low temperature. 
This is illustrated by Fig.^(a) where, for both densities studied, the solid line starts to deviate from the Monte Carlo 
data around f3 < 5. As in the repulsive case, this occurs because the self-consistent expression for double-occupancy 
Eq.(^) fails once we enter the renormalized-classical regime where a pseudogap appears. Following Ref. 1551 , we expect 
that the pseudogap in A( k f, uj) opens up as a precursor of Bogoliubov quasiparticles in the BSC ground state. Since 
this ground state starts to control the Physics, it is natural to expect that a high-energy quantity such as D should, 
in the pseudogap region, take the zero temperature BCS value, namely 

("t-i)-(^E-!) 2 + (^E^) 2 ' ( 9 ) 

k k 
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FIG. 2. (a) Double occupancy D calculated from the ansatz (solid curves), from the BCS ground state (dashed lines), and 
from the QMC simulations (circles) for \U\ = 4 at n = 0.5 and 0.8. (b) Renormalized interaction strength \U PP \ as a function 
of bare \U\ for T = 0.222 at n = 0.5. The dashed line is U pp = U. 

TABLE I. Double occupancy for various fillings n and temperature j3 obtained from QMC calculations and from the TPSC 
approach. The systematic error in the QMC calculation is of order 1%, larger than the statistical uncertainty 



U 



-4 

-4 



n 

0.8 
0.8 
0.55 
0.55 



4 
6 
4 
6 



QMC 
0.2671 ± 10~ 4 
0.2703 ± 10~ 4 
0.1561 ± 10~ 4 
0.1590 ±10~ 4 



TPSC 

0.273 

0.288 

0.157 

0.169 



where 
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with A the BCS mean-field gap. The chemical potential fi and gap A are determined self-consistently through 
the number and gap equations for given U , T and n. The value of D obtained with this approach is plotted as 
a dotted line in FigJ|(a), where it is apparent that the agreement with Monte Carlo calculations is excellent. In 
fact, for \U\ > 3, the agreement is always at the few percent level. For smaller \U\, deviations occur, probably 
because at small coupling the order parameter at q=0 does not dominate anymore the sum over all wave vectors in 
l/NY2$(A-(q)' A(—q)) = (A^A) = (n-j-ni). In fact, as we shall see below, a good estimate of double-occupancy may 
also be obtained at small U just from second-order perturbation theory. Now consider the temperature dependence 
of double-occupancy. The dependence predicted by the finite-temperature BCS result is on a scale T = 1 at \U\ = 4. 
That temperature dependence is clearly wrong for our problem since above Tmf the BCS approach would give us 
back the non-interacting value. Hence the BCS result may be used only in the following way. For \U\ > 3 we can 
use the T = value of {n^n^) in the pseudogap region and the self-consistent value Eq.(||) above it. In the general 
case, (ni-ni) does depend on temperature in the pseudogap regime, but that dependence should be relatively weak, 
as discussed in the repulsive case. Q 

We stress, however, that deep in the pseudogap regime our approach becomes eventually inaccurate when we obtain 
(n^ni) from the self-consistent equation (Q). The reason for the loss of accuracy is analogous to that found at n = 1 
in the repulsive |j2| case: In the present U < case, (n^ni) — > (nj,) as T — > to prevent a finite temperature phase 
transition. The approach also eventually becomes less accurate in the pseudogap regime when we take {n^n^} from 
BCS, but the fact that a more physically reasonable value of (n^ni) may be obtained in that case at T — helps 
extrapolate a little bit deeper in the pseudogap regime. The internal accuracy check discussed at the end of this 
section helps quantify the region of validity of the approach. 

Before moving on with the comparisons, a few comments on the actual renormalized interaction strength \U PP \ 
resulting from the Two-Particle Self-Consistent calculation. In Fig^fb) | U pp \ (denoted as stars) is plotted for T — 0.222 
as a function of bare U by using the self-consistent expressions Eqs.(||) and ([?]). For this temperature, U pp approaches 
bare [/for \U\ < 1 while in the intermediate coupling regime 1 < \U\ < bandwidth one notices a strong deviation 
from the bare \U\ . This deviation of U pp from U makes a drastic difference, in particular in the two-particle function 
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FIG. 3. Calculated s-wave paring structure factor S(q, r = 0) (filled triangles) and QMC S(q, r = 0) (circles) for \U\ = 4 and 
various temperatures, (a) at n = 0.5 and (b) at n = 0.8 on a 8 x 8 lattice. The dashed lines are to guide the eye.. 



Xp (s) which ultimately governs the particle dynamics, including the pseudogap behavior. The saturation of \U PP \ in 
Fig.g(b) signals the onset of the strong-coupling regime inaccessible in our approach that does not take into account 
the strong frequency dependence of vertices and self-energy present in this limit. 

Our approach gives a very accurate result for double-occupancy, but what about correlation functions? Another 
two-particle quantity related to the susceptibility is the pair structure factor. In Fig.H the calculated s-wave pair 
structure factor (filled triangles) S(q,r = 0) = (A(q)A^(-q) + A($+ A(-q)), where A{q)+ = -^Er c t t+ c Z > is 
compared with our QMC results (circles) for \U\ = 4 at n = 0.5 (top panel) and 0.8 (bottom panel) on a 8 x 8 
lattice. As the temperature decreases, the q = mode becomes more singular in both results, a characteristic feature 
for growing s-wave pairing fluctuations. In most of the Brillouin zone, the agreement is excellent, in particular, for 
n — 0.5 where the maximum difference is less than 10% (Fig. |](a)). For n — 0.8 our calculated structure factor 
overestimates QMC results at most by 20%. 

At the first level of approximation, we can also estimate the interaction-induced shift in chemical potential by 
starting from our approximate expression for the self-energy YS l \ Eq.(|5|). Let us call the corresponding chemical 
potential iJ- 1 ' = A*o + Tf- 1 ' . Our best estimate of single-particle quantities is obtained from the self-energy at the 
second level of approximation, Eq.(|[). The corresponding chemical potential /i^ 2 -* is calculated by requiring that the 
filling be the same as the one used at the first level of approximation. This procedure is identical to that for the 
repulsive model 56 and was suggested by Luttinger. It is discussed also in Sec.IV-C of the accompanying paper. |37]] 
Figure || illustrates how the two estimates for the chemical potential on a 8 x 8 lattice converge towards the value 
obtained from QMC calculations for the same size lattice. Consider the data for filling n = 0.6, in the upper part of 
the figure. The open squares, representing // 2 ) — /i , agree, within the error bars, with the QMC data represented by 
open circles. The first estimate for the chemical potential shift /i^ 1 ) — /iq (open diamonds) starts to deviate from both 
the QMC data (open circles) and from /i^ 2 ) — (1q (open squares) at the temperature where the pseudogap opens up 
(see Figjl^ below). Below this temperature, the self-energy becomes strongly frequency and momentum dependent, 
a feature captured by our improved estimate E 1 - 2 -' for the self-energy, but not by our first estimate, E' 1 '. The filled 
squares (/i' 2 ' — fJ-o) and filled diamonds (fi^ — fJ-o) were obtained for a 64 x 64 lattice. They illustrate that one should 
compare finite-size QMC calculations to many-body calculations done on same size systems. They also illustrate that 
the second estimate for the chemical potential shift (squares) is more sensitive to system size. This is expected from 
the fact that it is only at the second level of approximation that the pseudogap appears for large correlation lengths. 
Let us now move to the lower part of Fig.Q, where results closer to half-filling are plotted. The chemical potential 
there is very close to its exact temperature-independent half-filling value, U/2, hence there is not much room to see 
the difference between the first and second estimate for the chemical potential. Nevertheless, even for n = 0.95, the 
second estimate (open squares) is closer to the QMC data (open circles) than the first estimate (open diamonds). 

We now compare the momentum distribution, a static quantity, obtained from QMC and from our analytical 
approach at the second level of approximation S^ 2 ^. The momentum dependent occupation number n(k) (solid 
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FIG. 4. Comparison of the chemical potential shifts fj,^ 1 ' — fio (open diamonds) and /i' 2 ' — /in (open squares) with the results 
of QMC calculations (open circles), all done on a 8 x 8 lattice. QMC error is a few percent, smaller than the open circles. 
Upper set of points are for n = 0.6, and lower set of points are for n = 0.95. For n = 0.6, many-body calculations for a 64 x 64 
lattice are also illustrated by filled squares for — fxo , and by filled diamonds for fj,^ 1 ' — (jlo- In all cases, fi^ — fj,o is closer 
to QMC data than /i' 1 ' — fiQ. 
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FIG. 5. The momentum dependent occupation number n(k) for \U\ = 4 and n = 0.5 at T = 0.25. The circles denote n(k) 
from QMC calculations by Trivedi et al. |57]] on a 16 x 16 lattice. The solid curve is calculated according to the equations given 
in this paper, while the dashed one is computed by replacing U pp by U in the self-energy with all the rest unchanged. The 
long-dash line is the result of a self-consistent T-matrix calculation, and the dot-dash line the result of second-order perturbation 
theory. 
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FIG. 6. The single-particle density of states for \U\ = 4 and n = 0.87 at different temperatures on a 8 x 8 lattice, (a) T = 1/2, 
(b) T = 1/4, (c) T = 1/6, and (d) T = 1/8. Both the solid line obtained from many-body calculations and the dashed line 
taken from the QMC calculations of Ref. Jk| are obtained by analytic continuation of imaginary-time data using maximum 
entropy. The absolute error chosen for the Maximum Entropy continuation of the imaginary time many-body Green function 
for panels (a) to (d) is, respectively, 0.003,0.004,0.003,0.001. 



curve) is plotted in Fig. |B| along with the QMC calculations (circles) by Trivedi et al. |57j] for a 16 x 16 lattice, in 
a regime where the size dependence is negligible. The momentum distribution n(k) drops rapidly near the Fermi 
surface, corroborating that for \U\ — 4 the electrons are in the degenerate state, instead of in the non-degenerate 
state of the strong coupling regime or of the preformed-pair scenario. The agreement between the non-perturbative 
method and QMC is clearly excellent. If we had assumed a Migdal theorem and taken U 2 instead of UU pp in Eq.(||) 
for the self-energy, then we would have obtained the dotted curve, which in absolute value differs as much from the 
QMC result as a non-interacting Fermi distribution would. Clearly Migdal's theorem does not apply for this problem. 
In addition, Fig.|| also shows (long dashes) that a self-consistent T— matrix calculation does not compare to Monte 
Carlo data as well as our approach. Similarly, second-order perturbation theory (dot-dash) does not do well. 

In Fig. H we present the total density of states TV -1 Yl~~j? A(k ,u) for [7 = 4 and n — 0.87 on an 8 x 8 lattice and 

compare with existing QMC calculations by Moreo et al. [p8[ The many-body non-perturbative calculation is done in 
Matsubara frequency and analytically continued to real frequency using the same Maximum Entropy (ME) technique 
that is used for QMC calculations. [[59] This allows us to smooth out the results in the same way as in the QMC 
calculation, namely by including statistical uncertainties in our imaginary-time data. This procedure was discussed 
in Ref. |3^]. The uncertainties in the QMC data |58j that are presented in Fig.^ were not quoted. Hence, we chose the 
statistical uncertainties in the corresponding many-body calculations in such a way that the calculations presented in 
Fig.^(d) have the same degree of smoothness as the corresponding QMC data. More specifically, these uncertainties 
are of order 0.003 on the absolute value of the imaginary time Green functions, which is typical of Monte Carlo 
calculations. At T = 1/2 the density of states is similar to that for the non-interacting system. At T < 1/4, however, 
the spectral weight near the Fermi energy begins to be suppressed significantly with decreasing temperature, leading 
to a pseudogap. The small shoulders in the intermediate frequency regime for T = 1/6 and 1/8 come from finite size 
effect. When we use a 64 x 64 lattice, these shoulders completely disappear. 

In the case of the single-particle spectral weight, the latest QMC calculations at n — 1 |33| included studies of the 
finite-size effects, of the imaginary-time discretisation and of the uncertainties induced by the size of the Monte Carlo 
sample. They have shown that, at half-filling, there is indeed a pseudogap, in contrast to earlier findings ]36|j . Detailed 
comparisons with the many-body approach analog to the present one have been done. Although these studies were for 
the repulsive model, at half-filling the results apply for the attractive model since they are canonically equivalent at 
n = 1. One only needs to generalize the many-body approach to include the presence of the SO (3) symmetry, as done 
in Ref. J33|. Slightly away from half-filling, namely for n = 0.95, QMC simulations |5l| [[30| have found that a similar 
pseudogap opens up for U = —4 at a temperature large enough that the SO (3) symmetry present at half-filling is 
barely broken (Tx S> /«) . 

Here we present results for A( k , u>) and for the total density of states N (u>) = N^ 1 Y^T? A( k , u>) at n = 0.8, 
U = —4 and two different temperatures, T — 0.25 and T = 0.2, in the pseudogap regime. For that filling, finite-size 
studies are complicated by the fact that the available wave vectors do not necessarily coincide with the non-interacting 
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FIG. 7. Total density of states (first column) and single-particle spectral weight at wave vector (0, 37r/4) (second column) are 
shown for U — —4, n — 0.8 and two different temperatures. The lattice size is 8 x 8. Dashed lines are the result of Maximum 
Entropy continuation of QMC data, and solid lines the result of Maximum Entropy continuation of many-body results with the 
same errors added in. Plots (a) and (b) on the first line are, respectively, N(u>) the total density of states, and ^4(0, 37r/4, u) 
the single-particle spectral weight at a point near the Fermi surface, for T = 0.25. Panels (c) and (d) on the last line are the 
same plots but for T — 0.2. The number of measurements done is 1.2 x 10 5 for T = 0.25 and 1.6 x 10 5 for T = 0.2. In all cases, 
the absolute statistical error on the imaginary-time data is of order 2.0 x 10~ 3 . Singular values less than 10 -3 are dropped in 
the Maximum Entropy inversion. 



Fermi surface. The closest such wave vector, k — (0, 37r/4), is the one chosen for comparisons of A( k , ui). The results 
of QMC calculations are shown in Fig.(fjj) as dashed lines and those of the many-body approach as solid lines for the 
same 8x8 system size. The imaginary-time data corresponding to Fig. (Q) appears in Fig.([|). It is clear that all the 
main features of the Monte Carlo are seen by the many-body approach. The relative size of the split peaks at the 
Fermi surface is very sensitive to the actual location of the Fermi wave vector. 

We conclude this section with the accuracy check described in the accompanying paper |37j and in previous work 
p2| . The question here is whether it is possible to find the domain of validity of the approach even in the absence of 
Monte Carlo data. The answer is given by Table II, where all results were computed for a 64 x 64 lattice. One should 
focus on the three columns labeled Many-Body. In the first of these, one finds the value of U (rtjftj) computed with 
TPSC, i.e. from Eq.(|) and the local -pair sum rule Eq.(^|). That number is the same as that which would be obtained 
from Tr[£( 2 )G( 1 - ) ] . However, finding by how much it differs from Tr[E^ 2 ^G' 2 '] , listed in the second column labeled 
Many-Body, gives an indication of how much the theory is internally consistent. The third column labeled Many- Body 
gives the absolute value of the relative difference between the first two Many-Body columns and it helps appreciate 
where the theory fails. Clearly, as temperature decreases and one enters the renormalized classical regime, the theory 
becomes invalid since the difference between U {n^n^) and TrfE^G^ 2 -'] starts to increase rapidly. As expected also, 
the theory is better at smaller coupling. Instead of computing U {n^n{) with the TPSC approach Eq.(^), one may 
also take it from the zero-temperature BCS value as described above. Comparing with the corresponding Many-Body 
columns, it is clear the BCS estimate of double occupancy does not compare as well with the TPSC result at small 
coupling, as explained above. In fact, in this region, second-order perturbation theory compares better. This is partly 
because the Hartree-Fock contribution to double-occupancy is dominant. The worse results for perturbation theory 
are for quarter filling at U — —4. Near half-filling it is known that second-order perturbation theory works well again. 
Note also that even if we start from the BCS result for U (n^ni), one can tell that the theory is becoming less accurate 
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FIG. 8. Imaginary-time data from which are extracted the density of states and spectral weights appearing in the respective 
panels of Fig.f? . Error bars are shown only for the Monte Carlo data but the analytical data has been added the same errors 
as the Monte Carlo data. 



at too low temperature because, there again, U (n-\n\) and Tr [S^G^ 2 -*] begin to differ more and more from each 
other. This internal accuracy check, however, cannot tell us which approach is more accurate compared with QMC. 



III. PSEUDOGAP FORMATION IN THE DENSITY OF STATES AND IN A(K F ,u) 



In the hrst subsection below, we summarize pre vious work on pseudogap in the attractive Hubbard model. In the 
second subsection, we use the approach of section [I A to study the conditions under which a pseudogap appears in 
the density of states and in A{ k p, u>). 



A. Overview of some recent work 



The effect of superconducting fluctuations on the density of states was studied long ago. j61j To elucidate further the 
Physics of pseudogap formation, especially in the single particle spectral weight, many theoretical studies have focused 
on the attractive Hubbard model. An exhaustive review may be found in Ref. |Q. Although the attractive Hubbard 
model is clearly not a realistic model for cuprates since it predicts an s— wave instead of a d— wave superconducting 
ground state, it is an extremely useful paradigm. Indeed, except for the lack of cutoff in the interaction, it is analogous 
to the BCS model and is the simplest many-body Hamiltonian that leads to superconductivity with a possible crossover 
from the BCS limit at weak coupling to the Bose-Einstein limit at strong coupling. (6^] A key point, as far as we are 
concerned, is that it also represents the only Hamiltonian for which QMC simulations are available now as a means 
of checking the accuracy of approximate many-body calculations. 

The conventional picture [ p2| , based on mean-field ideas, is that in the weak coupling regime (|£7/t| <C 1) pairing 
and phase coherence happen at the same temperature while in the strong coupling regime (\U/t\ 1) phase coherence 
may occur at T c much lower than T* where pair formation happens. In the latter case, the Fermi surface is destroyed 
well before the superconducting transition occurs. Since the ARPES experiments suggest a relatively well-defined 
Fermi surface, it has been suggested |i3|,[57],[34| that at intermediate coupling it is possible to retain aspects of both 
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TABLE II. Internal consistency check, for various U, n, and T on a 64 x 64 lattice. The three columns labeled Many-Body 
give, respectively, the value of U < n^ni > obtained at the first level of approximation with the TPSC approach, that obtained 
at the second level of approximation from Trp' 2 ' G^] and the absolute value of the difference between the two results. The 
following column is Tr[EG] obtained from second-order perturbation theory, and the last three columns are analogous to those 
where double-occupancy is obtained from the TPSC but they start from the T = BCS estimate of U < n-|-nj >. 
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the Fermi surface of weak coupling and the preformed pair ideas of strong coupling. These possibilities have been 
extensively studied by several groups using a number of approaches that we will crudely divide in two types, numerical 
and many-body approaches. 

Previous numerical QMC work charted the phase diagram of the attractive Hubbard model fu1|f48| -[50| . They have 
also investigated the pseudogap phenomenon from intermediate to strong coupling ](36|-|68|] . We stress that the Physics 
in strong coupling is different from the weak to intermediate coupling limit we will study below. On the weak-coupling 
side of the BCS to Bose-Einstein crossover, there have been numerical studies of BKT superconductivity |58|,|30| as 
well as several discussions of pseudogap phenomena in the spin properties (susceptibility and NMR relaxation rate) 
and in the total density of states at the Fermi level. |33|,|7],[36| We have studied, through Quantum Monte Carlo 
(QMC) simulations, the formation of a pseudogap in A (kp,u>) in the weak to intermediate coupling regime of interest 
here p9|pl[ |. The present work is in agreement with our earlier results, as will be discussed in the next subsection. 

The many-body techniques that have been applied to the attractive Hubbard model in the weak to intermediate 
coupling regime are mostly T-matrix and self-consistent (Fluctuation Exchange Approximation) T-matrix approaches 
p9| fT2|. Let us consider the pseudogap problem in the non-superconducting state. At low density [ |73|j74| ], or with 
additional approximations j75|,[76|], a pseudogap may be found. By contrast, when the q — superconducting mode 
is relaxational, self-consistent T-matrix calculations J77J have failed to show a pseudogap in the one-particle spectral 
function A (kp, oj) = —2 Im G R (kp, us), a dimension-independent result. In two-dimensions, this absence of pseudogap 
is in sharp contrast with various QMC results ]39| , pT[ | and with general physical arguments inspired by studies of the 
repulsive case |55| [Q, which have already stressed that space dimension is crucial in the Physics of pseudogap 
formation. 

In the non self-consistent version of the T-matrix approximation however, a pseudogap can be found ]79| . Neverthe- 
less, since the T-matrix approximation takes into account the Gaussian (first nontrivial) fluctuations with respect to 
the "mean field state" , one important pathology of the T-matrix approximation in two dimensions is that the Thouless 
criterion for the superconducting instability occurs at a finite mean-field temperature Tmf- In the very weak coupling 
regime, this is considered inconsequential since the relative difference between Tmf and the Berezinskii-Kosterlitz- 
Thouless |8(| temperature (BKT) Tbkt is of order J8lj Tmf/Ef- In the intermediate-coupling regime however, this 
argument fails and one of the questions that should be answered is precisely the size of the fluctuating region where 
a pseudogap is likely to occur. 

By contrast with non self-consistent T— matrix calculations, self-consistent T-matrix approaches do have a large 
fluctuation region, but they assume a Migdal theorem which means they do not take vertex corrections into account in 
the self-energy formula and use self-consistent Green functions. In the same way as in the case of repulsive interactions 
fj2f fi3[ , the failure to treat vertex corrections and fluctuations at the same level of approximation may lead to incorrect 
conclusions concerning pseudogaps in A(k p,u). Self-consistent calculations have in fact lead to the claim |j35| that 
only d— wave superconductivity may have precursor effects in A {kp, lo) above the transition temperature while Monte 
Carlo calculations |59| |5l[] have exhibited this pseudogap even in the s— wave case. Recently, it has been pointed 
out using a different approach |86| that in the intermediate-coupling regime and when the filling is low, it becomes 
possible to have a bound q — pair. This Physics leads to a pseudogap in A(k, w) but it requires strong particle-hole 
symmetry breaking and is not specific to two dimensions. This result is discussed further in the following subsection. 

Some authors p9|- |2l| , p^ | have included phenomenologically a BKT fluctuation region either in T— matrix like calcula- 
tions, through Hubbard-Stratonovich transformation or otherwise [p2| . These calculations allow for phase fluctuations 
in the presence of a non-zero expectation value for the magnitude of the order parameter. In such a case, there is 
generally a real gap in A( k p, to), and additional effects must be included to fill-in the gap to transform it into a 
pseudogap |38[ . In the approach that we take, any SO (n > 2) theory would give qualitatively the same result above 
either Tbkt for n = 2 or above T = for n > 2. In addition, in our approach the magnitude of the order parameter 
still fluctuates. It is when both amplitude and phase fluctuations enter the renormalized-classical regime, i.e. become 
quasistatic, that a pseudogap may open up. 



B. Pseudogap formation in weak to intermediate coupling 



In this section, we focus on the Physics of fluctuation-induced pseudogap in the single-particle spectral weight 
of the two-dimensional attractive Hubbard model. This Physics has been discussed in prev ious QMC J39| , ^l| , |60[ 
and analytical work [ [55f but the present quantitative approach, based on the equations of Sec. II A, allows us to do 
calculations that are essentially in the thermodynamic limit and that can be done sufficiently rapidly to allow us to 
address other questions such as the crossover diagram in the temperature-filling plane. 

Since the cuprates are strongly anisotropic and may be considered as a quasi 2D systems, it is important to under- 
stand in detail the limiting case of two-dimensions. Mean-field theory leads to finite-temperature phase transitions 
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FIG. 9. (a) The density of states and (b) the spectral function at the Fermi surface for \U\ = 4 and n = 0.5 at different 
temperatures. The solid, dotted, dashed, long-dashed, and dot-dashed curves correspond to T = 0.143, 0.182, 0.211, 0.308, and 
0.364, respectively, except for panel (b) where the long-dashed line was calculated with a 2560 x 2560 lattice to illustrate that 
size effects are small. 



even in low dimensional systems where breaking of continuous symmetries is strictly forbidden at finite tempera- 
tures (Mermin- Wagner theorem |87|). Thus in low dimensions, mean- field theory, or fluctuation theory based on the 
mean-field state, lead to qualitatively wrong results at finite temperatures. One of the particular features of 2D that 
is captured by our approach, as we will see below, is that the mean-field transition temperature is replaced by a 
crossover temperature below which the characteristic energy of fluctuations is less than temperature, the so-called 
renormalized-classical regime. In this regime, the correlation length (£) increases exponentially until, in the super- 
conducting case, one encounters the BKT |8(| topological phase transition. As a result we find that when £ 3> Vp/T, 
the electronic system simulates the broken-symmetry ground state (£ = oo) at temperatures that are low but not 
necessarily very close to the transition temperature, leading to precursors of Bogoliubov quasiparticles ]5j| above 
T c . Recent dimensional-crossover studies using an analogous approach |7§| ] have suggested how the pseudogap will 
disappear when coupling to the third dimension is increased. 

The results of this section for two-particle properties, such as the pairing susceptibility and the characteristic pairing 
fluctuation scale £, are all computed for an effective 2560 x 2560 lattice. The one-particle properties are calculated 
instead for a 64 x 64 lattice in momentum space. Fast Fourier transforms (FFT) were used in that case to speed up 
the calculations. The solid and long-dashed lines in Fig.||(b), obtained respectively for a 64 x 64 and a 2560 x 2560 
lattice, illustrate that for single-particle properties a 64 x 64 lattice suffices. This lattice size, for the temperatures we 
consider, is large enough |3^] compared ^th (panel in Fig.[ll](b}) that single-particle properties are essentially the same 
as they would be in the infinite-size limit. Refs. |3lf| and p3| discuss how single-particle properties become rather 
insensitive to system size even when £ > L, as long as the condition ^ t h < L < £ is satisfied. This is discussed further 
in In the calculations, equation (Q) is solved iteratively, then the self-energy Eq.(||) is obtained in Matsubara 

frequencies. The analytic continuation from Matsubara to real frequencies are performed via Pade approximants fl90| . 
In order to detect any spurious features associated with this numerical analytical-continuation, we also performed 
real- frequency calculations. The results are identical, except for the fact that the Pade technique smooths out some of 
the spiky features of the real frequency formulation that are remnants of finite-size effects when the small imaginary 
part r\ in retarded propagators is very small. 

In Fig. ^ we show, for various temperatures, the total density of states (Fig. ||(a)) as well as the spectral function 
(Fig. ||(b)) A( k p, uj) for the Fermi surface point crossing the (0, 0) — (ir, 0) line for U = — 4 and quarter filling n = 0.5 
For T = 0.364 (dot-dashed curve) the density of states, on the left panel, is similar to that for non-interacting electrons. 
With decreasing temperature below T = 0.32, the low frequency spectral weight begins to be suppressed, leading to 
a pseudogap in the density of states. The condition for the appearance of a pseudogap in the spectral function 
A( k f,oj), on the right panel, is more stringent than that in the total density of states. Although the pseudogap in 
the density of states is well developed for T — 0.211 (dashed curve), it disappears in the spectral function for the 
same temperature. It is easier to form a pseudogap in the total density of states because of its cumulative nature: 
It suffices that scattering become stronger at the Fermi wave vector than at other wave vectors to push weight away 
from uj = 0. Hence, a pseudogap may occur in the density of states even if A( k p, u>) remains maximum at u) = 0. This 
is what occurs in FLEX (Fluctuation Exchange) type calculations. |6^j7^] It is more difficult to create a pseudogap in 
A( k f,u>) itself since, at this wave vector, transforming a maximum at uj — to a minimum requires the imaginary 
part of the self-energy to grow very rapidly as T decreases. J32] The generality of these arguments suggests that 
d— wave pairing fluctuations, which were considered in Ref. |{7Sj , ft8[ for example, should also lead to a pseudogap in the 
density of states before a pseudogap in A{ k f,uS). This feature is consistent with the recent experimental observations 
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FIG. 10. (a) The imaginary part of the pairing susceptibility at q — for T = 0.19 and (b) the characteristic low energy 
scale (stars) for pairing fluctuations at different temperatures for \U\ — 4 and n = 0.5. The inset in (a) is the imaginary part 
of the pairing susceptibility divided by frequency at q = for T = 0.19. 



H on high-temperature superconductors where pseudogap phenomena appear at higher temperatures in tunneling 
experiments than in ARPES experiments. Note also that with increasing temperature the pseudogap in both the 
density of states and the spectral function appears to fill instead of closing. This behavior is also in qualitative 
agreement with tunneling |l and with ARPES experiments HQ. All the above results are consistent with Monte 
Carlo simulations . In addition to having found a pseudogap in the density of states Q , Fig. (||) , the more recent 
Monte Carlo simulations done in the present and earlier papers |^] |[lj have also shown that a pseudogap may occur 
in A( k p, uj) even in s— wave superconductors, contrary to the claims of Ref. ]85| . 

Fig. |^(b) also shows one other qualitative result which is a clear signature of intermediate to strong-coupling 
systems, analogous to the signatures seen in optical spectra of high-temperature superconductors j83|. In changing 
T by about 0.03, from 0.21 to 0.18, the spectral weight rearranges over a frequency scale of order one, i.e. over a 
frequency scale about 30 times larger than the temperature change, and 5 times larger than the absolute temperature 
0.2. In weak-coupling BCS theory, by contrast, spectral weight rearranges over a frequency scale of the order of the 
temperature change. The frequency range for spectral rearrangement observed in Fig. |^(b) would be even larger if 
the coupling was stronger. ]84| ] This is a consequence of the fact that wave vector can be a very bad quantum number 
for correlated systems so that a momentum eigenstate can project on essentially all the true eigenstates of the system. 
The loss of meaning of momentum as a good quantum number and the corresponding spectral weight rearrangement 
over a large frequency scale happens suddenly in temperature in Fig. ^(b) because the correlation length becomes 
large at a rather sharp threshold temperature where the system becomes renormalized classical, as we now discuss. 

Let us then demonstrate that the opening of the pseudogap in the single-particle density of states occurs when 
the pairing fluctuations enter the renormalized-classical regime. [|9| In Fig. |l^(a) the imaginary part of the paring 
susceptibility at q = for T = 0.19 and the characteristic frequency v c for pairing fluctuations are shown for U = —4 
and n = 0.5. Since for the parameters studied here the q = mode is deep in the particle-particle scattering continuum, 
it has the characteristic frequency dependence of a relaxational mode, 1/(1 — iv/v c ) 1 that leads to a maximum in 
the imaginary part at some characteristic frequency v c . Even though we do not have perfect particle- hole symmetry, 
the Fermi energy is still large enough compared with temperature that Imx pp (0,^) jv is very nearly even (Inset in 
Fig. [l(](a)). For other temperatures the behavior is similar. In Fig. |l^(b) v c is plotted as a function of temperature. 
At high temperatures v c is larger than T but below T ~ 0.31 — 0.32 the characteristic frequency v c becomes smaller 
than T, signaling that we are entering the renormalized-classical regime. This phenomenon was also observed in 
QMC calculations. p9[ | In this regime, the thermal occupation number for pairing fluctuations is larger than unity. 
Clearly, the appearance of a pseudogap in the density of states in Fig. ^(a) follows very closely the entrance in the 
renormalized classical regime. 

The present results should be contrasted with those of Levin's group [pw . The pseudogap in their work comes from 
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FIG. 11. Spectral function for \U\ = 4 and n — 0.5 (a) along (0,0) — (0, tt) direction and (b) along (0,0) — (tt, tt) direction. 
The inset in (b) shows the pairing correlation length (solid curve), and £tu = vf/T along (0, 0) — (tt, tt) direction (dotted curve) 
and (0, 0) — (0, tt) direction (dashed curve) for the same parameters. 



the presence of a q = resonant pair state in the T-matrix. As the interaction strength decreases or the particle 
density increases, the q = bound state enters into the particle-particle continuum, thereby acquiring a finite lifetime. 
As long as the q = pair state is near the bottom of the scattering continuum it can remain a resonant state with 
a relatively long lifetime. Thus the origin of a pseudogap in their study is analogous to the preformed-pair scenario 
where the q = pair is separated from the scattering continuum. Such a resonance corresponds to strong particle-hole 
asymmetry in the imaginary part of the pair susceptibility. In order to have such an asymmetry for moderate coupling 
strength, very small particle density is required in this approach. In our case, the pseudogap occurs even when the 
particle-hole symmetry is nearly perfect. Furthermore, in our case, other factors like density and interaction strength 
do not influence the results in any dramatic way. Low dimensionality is the key factor since phase space is behind the 
existence of both the renormalized classical regime and the very strong scattering of electrons on the corresponding 
fluctuations. The ratio £/£th controls the importance of this scattering |3l[] as we discuss in the following paragraph. 

In Fig. [ll] we contrast the onset of the pseudogap in the spectral function on the Fermi surface along different 
directions, namely the (0,0) — (0, tt) and (0,0) — (7r,7r) directions, for \U\ — 4 and n = 0.5. At this density, where 
the Fermi surface is nearly circular, the anisotropy happens in a very small temperature range around T = 0.19. For 
T = 0.19 (dotted curves), the figure shows that the pseudogap occurs only along the (0,0) — (0, tt) direction. This 
anisotropy of the pseudogap in the spectral function should be contrasted with the fact that in the superconducting 
state, the gap is isotropic. The anisotropy at the temperature where the pseudogap opens up can be understood 
following the arguments of Vilk tt al. |32). Using the dominant renormalized-classical fluctuations (iq n = 0), these 
authors showed that for £ 3> vp/T the scattering rate (imaginary part of the self-energy) on the Fermi surface becomes 
large, leading to a minimum in the spectral function at u> = instead of the maximum that exists in the absence of 
a pseudogap. In the inset of Fig.|ll](b) the pairing correlation length £, as well as £ t h — vf/T along the (0, 0) — (tt, tt) 
and (0, 0) — (0, tt) directions are plotted as a function of temperature. Clearly £ grows exponentially with decreasing 
temperature. Furthermore, according to the above criterion, a pseudogap in the spectral function exists along one 
direction and not along the other when £ (solid curve) is larger than along (0, 0) — (0, tt) (dashed curve) but smaller 
than £th along the (0,0) — (tt, tt) (dotted curve), namely, in the temperature range 0.175 < T < 0.185. We obtain 
quantitative agreement with Fig. |ll](b) if we use £ = l.3£, t h as the criterion for the appearance of a pseudogap. While 
the pseudogap anisotropy happens in a narrow temperature range at this density due to the small Fermi velocity 
anisotropy (about 1.35), closer to half-filling it occurs in a large temperature interval since the Fermi velocity is nearly 
vanishing close to the (0, tt) point. |nj 

Finally in Fig. [l^ we present the crossover diagram for the pseudogap in the 2D attractive Hubbard model for \U\ = 
4. The dotted curve is a rough QMC 0] estimate (probably an upper bound) for the BKT transition temperature 
Tbkt- For all densities a pseudogap in the one-particle functions appears in a wide temperature range Tbkt < T < T* 
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FIG. 12. The crossover diagram of the 2D attractive Hubbard model for \U\ = 4. The filled triangles and stars denote the 
temperatures where a pseudogap appears in the density of states and the spectral function, respectively. The solid lines are 
a guide to the eye. The dashed curve is the BCS mean field temperature Tmf and the dotted curve is an estimate of the 
Kosterlitz-Thouless temperature Tkt extracted from QMC results by Moreo et al. mm. 



where T* is typically several times of Tbkt — T c . The pseudogap occurs earlier in the density of states than in the 
spectral functions for most of densities. Near half-filling, however, the pseudogap appears more or less at the same 
temperature in the density of states and the spectral functions. In QMC for small systems, there seems to be a 
difference in the temperatures at which the two pseudogaps open up. j30| Performing a calculation with finite second 
neighbor hopping t' , we have confirmed that this almost simultaneous opening of the pseudogaps happens because 
of the strong influence of the Van Hove singularity, which leads to vp = 0, and not because of nesting. Finally, 
note that at half-filling one has perfect O (3) symmetry in this model so that the transition temperature vanishes, 
as dictated by the Mermin- Wagner theorem, while the pseudogap temperature continues to be large, following the 
trend of the mean-field transition temperature instead of that of the Tbkt curve. This shows that symmetry of the 
order-parameter space contributes to enlarge the temperature range where the pseudogap occurs, as expected from 
the corresponding enlargement of the renormalized-classical regime. |5l| 



IV. CONCLUSION 



In weak to intermediate coupling, the attractive Hubbard model can be studied quantitativel y w ith a non- 
perturbative approach |37j that directly extends the corresponding method for the repulsive model ||4l|,pl|,p2[ , The 
simple equations of section II A are all that needs to be solved. This many-body approach has an internal accuracy 
check, no adjustable parameter and it satisfies several exact sum rules fl37|| . We have demonstrated the accuracy of this 
method through detailed comparisons of its predictions with Quantum Monte Carlo simulations of both single-particle 
and two-particle correlation functions. 

On the Physical side, we studied the fluctuation-induced pseudogap that appears in the single-particle spectral 
weight, in agreement with Monte Carlo simulations and in close analogy with the results found before in the repulsive 



case 1 32 [B3|. A key ingredient for this pseudogap is the low dimensionality. Indeed, in two-dimensions the finite- 
temperature mean-field transition temperature is replaced by a crossover to a renormalized-classical regime where the 
characteristic pairing frequency is smaller than temperature and the pairing correlation length £ grows faster than the 
single-particle thermal de Broglie wavelength ^ t h- In this approach, where vertex corrections and Green functions are 
taken at the same level of approximation in the self-energy expression |32]| [|33"f , the renormalized-classical fluctuations 
and the relatively large phase space available for them in two-dimensions lead to precursors of the superconducting 
gap (or Bogoliubov quasiparticles) in the normal state. This pseudogap can occur without resonance in the pair 
susceptibility pr| , and it appears not only in the total density of states but also in the single-particle spectral weight, 
in sharp contrast with what was found with self-consistent T— matrix approaches fn\ . Our approach fails at strong 
coupling or at low temperature very close to the Berezinskii-Kosterlitz-Thouless (BKT) transition. 

For \U\ — 4, the pseudogap regime occurs over a temperature scale that is several times the BKT transition 
temperature. The crossover to the renormalized-classical regime is about a factor of two lower than the mean-field 
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transition temperature but it has the same filling dependence, which can be quite different from that of the real 
transition temperature, which is strongly dependent on the symmetry of the order-parameter space |5lf . It is clear 
also that SO (2) (or U (1)) symmetry is not essential to the appearance of a pseudogap. It would also appear if there 
happens to be a hidden continuous symmetry group |92| |5l|] SO (n) with n > 2 describing the high-temperature 
superconductors. 

As stressed earlier in this paper, the attractive Hubbard model is not directly applicable to the cuprates. Neverthe- 
less, it helps understand the nature of superconducting-fluctuation induced pseudogaps, if they happen to be present. 
The pseudogap appearing for the underdoped compounds at high temperature in thermodynamic and transport mea- 
surements, or at high energy in tunnelling j(| and ARPES experiments, is most probably not of pure superconducting 
origin |93| l [Q. Nevertheless, close enough to the superconducting transition, in both the underdoped and over- 
doped regions, there should be an effective model with attraction describing the low energy Physics. Since even 
the high-temperature superconductors have a gap to Fermi energy ratio that is small, this effective model could be a 
weak-coupling one (but not necessarily [ p5| ) . Time-domain transmission spectroscopy experiments [ p6| in the 100GH z 
range suggest that the renormalizcd classical regime for the BKT transition has been observed in underdoped com- 
pounds, 10K to 15K above T c . Also, in the overdoped regime, recent experiments on the magnetic field dependence 
of NMR Tj -1 and Knight shift |97) suggest that the pseudogap appearing a few tens of degrees above T c is indeed 
a superconducting-fluctuation induced pseudogap. The pseudogap that we have described should appear in these 
regimes if an effective weak to intermediate-coupling attractive-interaction model is valid near T c . In this context, 
some of the important results that we found and explained are as follows. In the attractive Hubbard model the 
pseudogap appears earlier in the density of states than in the spectral function that would be measured by ARPES, 
as summarized in Fig.(|l2|). We also found, Fig.(^), that with increasing temperature, spectral weight appears to fill 
in the pseudogap instead of closing it. Finally, we also showed that as the system enters the renormalized-classical 
regime, spectral weight can rearrange over a frequency range much larger than the temperature scale. This is gener- 
ally a signature that momentum is becoming a very bad quantum number. Hence, for a given temperature scale, the 
frequency range over which the spectral weight can rearrange becomes larger with increasing coupling . All these 
features carry over the d— wave case |76l , |88fl . Qualitative differences between weak- and strong-coupling pseudogaps 
have been discussed in Ref . Q . 
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